Viable inflationary models ending with a first-order phase transition 
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We investigate the parameter space of two-field inflation models where inflation terminates via a 
first-order phase transition causing nucleation of bubbles. Such models experience a tension from 
the need to ensure nearly scale invariant density perturbations, while avoiding a near scale-invariant 
bubble size distribution which would conflict observations. We perform an exact analysis of the 
different regimes of the models, where the energy density of the inflaton field ranges from being 
negligible as compared to the vacuum energy to providing most of the energy for inflation. Despite 
recent microwave anisotropy results favouring a spectral index less than one, we find that there are 
still viable models that end with bubble production and can match all available observations. As a 
by-product of our analysis, we also provide an up-to-date assessment of the viable parameter space 
of Linde's original second-order hybrid model across its full parameter range. 

PACS numbers: 98.80.Cq 



I. INTRODUCTION 

One of the open questions in inflationary cosmology is 
the mechanism by which inflation came to an end. The 
current literature is dominated by two paradigms, viola- 
tion of slow roll bringing inflation to an end while the field 
is still evolving, and a second-order phase transition of 
hybrid inflation type. However, Guth's original (but un- 
successful) proposal [l[ invoked a first-order phase transi- 
tion whereby inflation ended by nucleation of bubbles of 
true vacuum. First-order transitions have subsequently 
experienced bursts of popularity. In the late 1980s, La 
and Steinhardt [2] initiated intensive investigation of 'ex- 
tended inflation' models, where modifications to Einstein 
gravity allowed bubble nucleation to complete in single- 
field inflation. A few years later those models were strug- 
gling in face of observations, and focus instead returned 
to Einstein gravity, now in a two-field context with one 
rolling and one tunnelling field 0, 0, HI, although see 
0.0- 

In addition to the usual quantum fluctuation mech- 
anism, first-order inflation models produce density per- 
turbations through the bubble collisions and subsequent 
thermalization. The spectrum of bubble sizes produced 
must be far from scale invariance to avoid clear violation 
with observed microwave anisotropies — the largest of 
the bubbles would otherwise be blatantly visible [sfilUl- 
This requirement is typically at odds with the need to 
maintain scale invariance in the spectrum produced by 
quantum fluctuations, a tension sufficient to exclude ex- 
tended inflation variants except in extremely contrived 
circumstances [ll|. The purpose of this paper is to in- 
vestigate whether the strengthened constraints of the 
post Wilkinson Microwave Anisotropy Probe (WMAP) 
era have eliminated the Einstein gravity first-order mod- 
els too and, by implication, assess whether it is plausible 
that voids exist below current detection limits. 

In Guth's original model, with one field, the inflaton 
must remain in the metastable vacuum long enough to 
allow for sufficient e-folds of inflation but in this case 



inflation never ends, the bubbles never thermalize and 
the transition doesn't complete. Introduction of a second 
field allows a time-dependent nucleation rate, permitting 
enough inflation to occur while the nucleation rate is low 
and a successful end when the rate rises to high enough 
values. This idea was proposed independently by Linde 
Q and, in more detail, by Adams and Freese [j under 
the name 'double-field inflation'. 



Typically the second field, which is trapped in the 
metastable vacuum, also provides most of the energy den- 
sity for inflation, although this depends on the particular 
values of parameters chosen. In that regime, the usual 
prediction is for a blue spectrum of density perturba- 
tions, ng > 1. In the last few years the trend in cosmic 
microwave background (CMB) observations has been a 
tightening of the confidence limits around a central value 
ns smaller than one, disfavouring this regime. Since our 
goal is to investigate the general viability of this type of 
model we will probe the entire parameter space, includ- 
ing the intermediate region where the contributions of 
each field to the energy density are comparable, making 
no approximations based on inflaton or the false vacuum 
domination. 



As stated above one expects these models to run into 
difficulty with recent observations closing in on a nearly 
scale invariant scalar spectrum. CMB anisotropies ob- 
servations place constraints on the maximum size of bub- 
bles that survive from a first-order phase transition, at 
the time when scales of cosmological interest leave the 
horizon. In turn this places a strong upper limit on the 
nucleation rate at this time, after which it must rise suffi- 
ciently to complete the transition and provide a graceful 
exit for inflation. In order to meet these two requirements 
the field must proceed swiftly along the potential, what, 
in light of observations, places the model under stress. 
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FIG. 1: (a) The potential for a second-order phase transition. The field reaches the true vacuum through a continuous 
transition, and the breaking of the symmetry implies that there will be defect formation at the end of the transition. The 
true vacuum minima develop once the field passes the point of instability, 0i ns t (b) The same for the first-order case. In this 
case the transition is discontinuous and proceeds through quantum tunneling of the ip field to the true vacuum. The second 
minimum develops after the point of inflection <^i n fl. The couplings in both (a) and (b) have been chosen so as to produce a 
visible barrier height (in working models this is negligible compared to the false vacuum energy). 



II. THE FIRST-ORDER MODEL 

We consider throughout a fairly general form of the 
potential for a first-order phase transition, given by 
Copeland et al || . 

V(cj>, ip) = ^A(M 4 + V 4 ) + ^aM 2 i/j 2 - \lM^ 

+^¥ + JaW. (i) 

This extends the simplest second-order hybrid inflation 
model by addition of the cubic term for the ip field. As in 
conventional hybrid inflation, one envisages that initially 
the inflaton field <p is displaced far from its minimum, 
and the auxiliary field ip is then held in a false vacuum 
state by its coupling to the inflaton. Perturbations are 
generated during this initial phase as cp rolls slowly along 
the flat direction. The dynamics in this region are pretty 
much those of single-field slow-roll inflation, though the 
auxiliary field ip may provide most of the energy density 
for inflation, see Fig. [1] 

In a model where the phase transition is second-order, 
shown in Fig.[T^,, the false vacuum becomes unstable after 
4> passes a certain value, </>; ns t, and the fields evolve clas- 
sically to their true vacuum (here producing topological 
defects as causally separated regions make independent 
choices as to which minimum to finish in). Although 
not the main topic of this paper, we explore current con- 
straints on this model in the Appendix. 

In the first-order case, shown in Fig. QJ>, if the pa- 
rameters in Eq. ([T]) are chosen appropriately, a second 



minimum develops once the field evolves past a point of 
inflection, <j>- m fi. At this point bubbles of the true vac- 
uum begin to nucleate and expand at the speed of light. 
The percolation rate is initially very small as the vac- 
uum energies are comparable, but as <fi approaches zero 
the interaction between the fields triggers a steep rise in 
the bubble production. Inflation ends when the nucle- 
ation rate reaches high enough values that the bubbles 
percolate and thcrmalize. In this case there is only one 
true vacuum and hence no topological defects. The chan- 
nel in which the field rolls after tunnelling is much too 
steep to sustain any inflation within the bubbles. 

For large values of <p there is only one minimum of 
the potential, and in the ip direction the potential looks 
like Fig. [2^,. However if j 2 > 4aA, a second minimum 
develops after (p reaches a point of inflection 



as in Fig. [2}d. The presence of the cubic term in the 
potential then breaks the degeneracy between the two 
minima, making it possible for the field to tunnel to the 
newly formed minimum. It is this second minimum that 
eventually becomes the true vacuum and the ip field be- 
gins to tunnel once the transition becomes energetically 
favourable, Fig. 

As mentioned in the previous section the quantum gen- 
eration of perturbations occurs away from this minimum, 
while the inflaton is rolling in the <fi direction, and we 
consider horizon exit to occur around 55 e-folds before 
the end of inflation [12j. This evolution of <p is a cru- 
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FIG. 2: (a) At early times, away from </> — 0, there is only minimum available for ij) and the field is trapped in the false 
vacuum, (b) Given appropriate choices for the couplings a second minimum begins to develop when the field reaches the point 
of inflection of the potential, (c) Once the transition becomes energetically favourable the ip field begins to tunnel to the newly 
formed minimum, which eventually becomes the true vacuum. 



cial feature of the model since it is the introduction of 
a time dependence in the tunneling rate that will allow 
the phase transition to complete, bringing inflation to an 
end. 

The rate at which the bubbles nucleate is given by the 
percolation parameter (the number of bubbles generated 
per unit time per unit volume), 



P = 



(3) 



In the limit of zero temperature (taken because the tran- 
sition occurs during inflation) the nucleation rate of bub- 
bles can be approximated by (l3| . 



P : 



XM 4 
4iJ 4 



exp(-5 E ) , 



(4) 



where Se is the four-dimensional Euclidean action. Se 
was obtained for first-order transition quartic potentials 
by Adams |14|, who fitted the result as 

4tt 2 

S E = —(2-5)- 3 ( ai S + a 2 6 2 + a 3 5 3 ), (5) 

where ai = 13.832, a 2 = -10.819, a 3 = 2.0765, and 6 
is a monotonic increasing function of </> 2 , 



9\a 9AA' 



7 2 7 2 M 2 



(G) 



The allowed range has < S < 2 (outside this range solu- 
tions correspond to energetically disallowed transitions) . 

The transition to the true vacuum is complete once the 
percolation reaches unity, (one bubble per Hubble time 
per Hubble volume), allowing the bubbles of the true 
vacuum to coalesce. 

However in the most general case inflation need not end 
through bubble nucleation. If the potential is too steep 
slow-roll is violated before bubbles thermalize and infla- 
tion ends before the transition completes. In this case the 
precise mechanism which completes the transition is ir- 
relevant given that it occurs after inflation ends, and for 



our purposes the scenario is indistinguishable from the 
single-field case. (In this paper we do not consider grav- 
itational waves produced via bubble collisions, but these 
may provide a further observable [ID, [H, 0, EH that can 
ultimately be used to constrain this type of model.) 

The distinction between the two possibilities is given 
by the two values of the field, that at which the nucleation 
rate reaches unity, and that which makes e ~ 1 (violation 
of slow-roll), where e is the usual slow roll parameter 
defined in Eq. @. Inflation ends by whichever value of 
<f> is reached first, 



<^end = max(^ e , cr it) ■ 



III. INFLATIONARY DYNAMICS 



(7) 



A. Regimes 

Two different regimes can be distinguished, regarding 
which field we wish to have dominate the energy den- 
sity. In the usual hybrid inflation regime the energy den- 
sity of the potential is dominated by the false vacuum 
AM 4 ^> m 2 </> 2 , which provides the energy for inflation. 
In the opposite regime, in which the inflaton dominates 
the energy density, the dynamics rapidly approach those 
of single-field inflation since, as we will see, slow-roll vi- 
olation occurs sooner. 

Working in either of these two regimes would allow us 
to simplify some of the expressions governing the dynam- 
ics during inflation, such as the number of e-folds and the 
slow- roll parameters, Eqs. (fT2"]l . © and (fT5")) , and to pro- 
ceed via an analytical treatment instead of a numerical 
one. However our purpose here is to probe the dynamics 
of the full ng — r parameter space, (r is the tensor-to- 
scalar ratio given by Eq. (fT3| ) so as to determine whether 
there still remain models consistent with CMB observa- 
tions. Hence we also include the intermediate regime in 
our analysis, where the energy densities of the two fields 
are comparable, particularly when the transition between 
slow roll violation and bubble nucleation occurs. For this 
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reason we will retain the full form of the potential and 
proceed through numerical calculations. 



B. Field dynamics 

In order to specify the dynamics of each model we be- 
gin by finding the field value, max , at which inflation 
ends so we need to determine (j) £ and C rit- 4>e is ob- 
tained by evaluating the first slow-roll parameter for our 
potential and taking it to unity, 1 



'"-pi 
16tt 



4 i ? 2 

m <p m 



pi 



tt(AM 4 + 2to 2 2 ) 2 



(8) 



Inverting for (j) yields, 



'pi 



± mmpiv/m 2 ™ 2 ! - 8ttAA/ 4 - AirXM 4 



87rm 2 

(9) 

and we take the largest value of <fi. Note that the solution 
exists only for large values of to, where m 2 mp[ > 87rAAf 4 . 

To determine cr it we need to find the value at which 
the percolation parameter reaches unity, p cr it ~ 1. Solv- 
ing Eq. ([3]), we get 



In 



AM 4 



4pcrit# 4 ' 



(10) 



where Sent is given by Eq. ([5]). 

Inverting Eq. (jTUJl yields a value for cr j t (only one of 
the three roots lies in the allowed range) and in turn this 
allows us to determine on( j, and, by comparison with e , 
the mechanism by which inflation ends. 

Knowing on( j we can calculate the value of the field 
at horizon exit, <p^. In this model <j) rolls towards its 
minimum at <p = so ^55 > 4> C nd- Using the expression 
for the number of e-folds between two field values <pi and 
4>2 we get, 

tf(fc,fc)=]n£~-i£. (11) 



CLl TO 



PI J<t>l 



For </>i = </>55 and 4>2 = 4>end, and substituting for V, we 
have 



N(4> 55 ,4> cn d) = 27TA- 



M 4 



TO 2 TO pl ln ^ + ^ 5 -^ d) ' 

(12) 

where we make no assumptions on the relative size of the 
two masses and retain both terms. Substitution of <p on d 
yields ^55 and now we can calculate the scalar spectral 



2tt 



1 The field -0 sits in the false vacuum during the inflationary phase, 
since this is the only minimum available to in this region of 
the potential. This happens regardless of the means to ending 
inflation, so ip is set to zero throughout this section. 



index, ns, and the tensor-to-scalar ratio, r, at horizon 
exit, by use of their expressions in terms of the usual 
slow-roll parameters, 



ns — 1 = — 6e + 2t] : 
r = 16e , 

where e is given by Eq. ([9]), and 77 is 



TOpi V" 



2 ™2 



TO TO 



PI 



8?r V 2tt(AA/ 4 + 2to 2 ^ 2 ) ' 



(13) 
(14) 



(15) 



where the last equality is obtained by substitution of the 
potential. 

At this point we can locate the model in the ns — r 
plane and determine its position in relation to WMAP5 
confidence limits 1191] . 



C. Choosing parameters 

Throughout we set the self-interaction and coupling 
constants, A and A' respectively, equal to unity. We are 
then left with two constants, a and 7, and requiring the 
energy density of the true vacuum to be zero fixes one of 
these in terms of the other. We will fix a in terms of 7 
but the reverse option could just as well be taken. 

The CMB amplitude normalization can be used to re- 
late the two masses. We use this to fix the mass of the 
light field <p and then we are left with only two undeter- 
mined parameters: the energy of the false vacuum, M , 
and the constant a. For each value of a, varying M fully 
determines the dynamics of the fields, and describes a 
trajectory in the ns — r plane shown in Fig. [3l 

Each line is composed of two branches which corre- 
spond to the two solutions of the WMAP normalization, 
and converge for large values of M ~ 2.7 x 10~ 3 mpi. For 
values of M larger than this there is no solution to the 
amplitude normalization hence no viable models. This 
can be seen also in Fig. 2] which illustrates how the two 
different approximation schemes converge to a common 
behaviour and cease to exist after a certain value of M 
(c.f. Fig. 1 of Ref. [|). 

The right-hand branch in Fig. [3] corresponds to the 
lower branch in Fig. [4] and to the smaller value of rn 
from the WMAP normalization. In this branch the ap- 
proximate relation M ~ to 2 / 5 (in Planck units) holds and 
the false vacuum dominates. The dynamics are indistin- 
guishable in the ns — r plane when M < 10 -4 mpi. We 
start with the typical slightly blue tilted spectrum and 
negligible tensor fraction. As to continues to increase so 
does the deviation from ns ~ 1 until the approximate re- 
lation between the two masses breaks down and we have 
the inflaton playing a more significant role in the relative 
contribution of the two fields. At this point we observe 
a turn in the ns — r plane, and the solution enters the 
intermediate region of comparable field energy densities. 

Despite this we still observe inflation ending by bubble 
nucleation throughout this branch, from small values of 
M to the maximum at M ~ 2.7 x 10 _3 TOpi. 
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FIG. 3: The trajectories described in the ns — r parameter plane for first-order models when M is varying and m is set by the 
CMB normalization. The two lines correspond to different values of the coupling constant a, the outermost a — 0.1 and the 
innermost a — 0.01. The two endpoints correspond to the endpoints in Fig. |4]and converge at M ~ 10 _3 mpi corresponding to 
the union of the branches in Fig. [4] at (ns, r) ~ (0.99, 0.5). Mass values are given in Planck units. 
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FIG. 4: The relation between the two mass scales. The 
WMAP normalization admits two solutions for m, corre- 
sponding to false vacuum domination over the inflaton (lower 
branch), and the opposite regime, for large m, which is nearly 
independent of M (upper branch). The two regimes converge 
to common behaviour. We consider all three regimes in the 
analysis and set a — 0.1, given by the upper curve in Fig. [3] 



In the opposite branch, on the left-hand side, the model 
starts inside the WMAP5 95% confidence contour, well 



inside the inflaton dominated regime. Similarly to the 
other branch we observe an initial period where there is 
little dependence on the false vacuum energy, correspond- 
ing to the plateau on Fig. 01 and the dynamics are very 
well approximated by those of standard single-field infla- 
tion with a <t> 2 potential, well known to satisfy WMAP5 
data. 



This regime breaks down as the false vacuum energy 
increases and eventually we recover the regime where 
the phase transition triggers the end of inflation be- 
fore the violation of slow roll, meaning we are again in 
the bubble production scenario. The interesting results 
here draw from the fact that the transition occurs inside 
the WMAP5 95% confidence contour, making these vi- 
able models even away from false vacuum domination. 
Fig. [5] is a zoom of this region showing the field mass, 
M, at which the transition to bubble nucleation occurs, 
A/ ~ 5 x 10 _4 mpi, still allowed by the 95% confidence 
limits. 
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FIG. 5: Zoom of Fig. [3]showing the transition between models ending by slow roll violation and bubble inflation. The transition 
happens around M ~ 5 x 10 -4 mpi, well within the WMAP5 95% confidence contour. Mass values in Planck mass units. 



IV. THREE CONSTRAINTS 

In the previous section we looked at constraints in the 
n S ~ r plane. By specifying a value for a, one of our 
two free parameters (M, a), the CMB normalization then 
allows us to recover a trajectory in this plane and as- 
sess where the density perturbations are compatible with 
WMAP5 data. We now compute other constraints on the 
scenario, in the M — a plane. 

A. Model consistency 

We begin with the requirement that M be not larger 
than an upper limit above which, for a particular choice 
of couplings, the transition does not complete (</> cr it does 
not exist). We call this the model consistency constraint, 
which translates to a relation for the value of (/> cr it , com- 
ing from the requirement that there exists a solution of 
Eq. ([6]) for S. Because of the constant term in Eq. ([6]) 
this is an additional requirement to < S < 2. 

Since we have chosen to set m by the CMB normal- 
ization this can be translated into an excluded region in 
the (M, a) plane (although alternatively we could have 
expressed it in terms of a region in (M, m) , by having a 



specified by the CMB normalization instead) . This yields 
the region below the upper (blue) curve in Fig. [(J] We see 
that specifying a value for the false vacuum density im- 
poses an upper limit on the coupling a (alternatively on 
the inflaton mass, m) in order for the model to have the 
possibility to complete the phase transition. 



B. Big bubble constraint 

We adopt here a fairly crude criterion to judge whether 
the bubbles are compatible with observations, which is 
that any bubbles produced at the end of inflation and 
expanded to astrophysical sizes must, during the epoch 
of recombination, have a comoving size not larger than 
20/i~ 1 Mpc 8]. This corresponds to a maximum filling 
fraction at that time of 10 -5 , and puts an upper bound 
on the percolation rate of bubbles at the time the scales 
we observe today left the horizon: 
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With our form for the action Eq. ([5]) and choice of po- 
tential this becomes 

S 55 ~ -2.9 + 4 In ^^ + 11.5. (17) 

This gives us the region between the short dashed (black) 
lines in Fig. [5] 

C. WMAP constraint 

We can similarly place constraints on the (M, a) plane, 
by considering the 95% confidence limit resulting from 
the WMAP5 ns — r plane when tensors are included. 

rc. s <1.05. (18) 

Inverting Eq. (|18[) gives us an upper limit on M in terms 
of a, resulting in the region left of the long dashed (red) 
line in Fig. [5] 

We also see from Fig.|6]that this constraint is opposed 
to that coming from the CMB maximum bubble size re- 
quirement, as we argued in Section [IJ Big bubbles at last 
scattering put an upper limit on the nucleation rate at 
horizon crossing while CMB constraints on the spectral 
tilt put a lower bound on the nucleation rate, from the 
requirement that ns is not too distant from scale invari- 
ance. 

Nevertheless, a region of parameter space survives all 
constraints. 



V. CONCLUSIONS 

Our principal conclusion is that there do remain Ein- 
stein gravity models of first-order inflation which are 
compatible with observations, despite the increasing ten- 
sion between the need for a scale-invariant primordial 
spectrum and the suppression of large-scale bubbles. We 
have exhibited a particular class of model and found the 
parameter region where the first-order model is viable. 
Its predictions for ns and r are similar to the simple 
m 2 (j) 2 slow-roll inflation model, though a little further 
from scale-invariance. 

In this paper we have imposed a relatively simple con- 
straint on the bubbles, and have then assumed that their 
impact on the CMB is negligible as far as constraints 
on the primordial perturbations are concerned. A more 
detailed treatment would combine the two perturbation 
sources and refit to the CMB data, which may lead to 
some modification to the outcome in regimes where the 
bubble production is close to the observational limit. For 
models where the bubbles are safely within the observa- 
tional limits this is not an issue. 

This paper demonstrates that we are still some way 
from having a clear view as to how the inflationary period 
of the Universe may have ended. The literature contains 
three different mechanisms — violation of slow-roll, a 



second-order instability during slow-roll, and bubble nu- 
cleation — and we have shown that the last (and least 
popular) of these remains a viable option. First-order 
models are of phcnomcnological interest as the bubble 
spectrum is an additional source of inhomogeneity that 
could be considered in matching high-precision observa- 
tions. The bubble collisions may also generate detectable 
gravitational waves [15j, |l6|, [171, |l8| . There is therefore an 
ongoing need to refine understanding of the nature of 
perturbations induced by a primordial bubble spectrum. 
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APPENDIX A: THE SECOND-ORDER MODEL 

Although not part of our main study, the full param- 
eter range of the second-order hybrid inflation model 
[E US HH is easily studied using the machinery we have 
used for the first-order case. The second-order model 
also uses the energy density of an auxiliary field to raise 
the energy scale for inflation without endangering slow 
roll. The phase transition in this case is continuous, 
with the ip field rolling down to the true vacuum (see 
Fig. HJa)). There are no bubbles now and hence no bub- 
ble constraint; we just have to consider whether the usual 
perturbations are compatible with WMAP5 data. Fur- 
thermore since now there is no cubic term to break the 
degeneracy between the two minima, there is the possibil- 
ity of topological defect formation at the end of inflation, 
as different regions in space roll towards one or the other 
minimum. However we do not consider their possible 
impact here. 

The dynamics are closely related to those in the first- 
order case. The critical point where the phase transition 
completes is a point of instability 0i ns t , after which %jj = 
becomes unstable and starts to roll. The potential for 
this case is a particularization of the first-order potential 
Eq. (fT]) with A = —a and 7 = 0, and becomes, 

U(0, = \\{^ 2 - M 2 ) 2 + + ±A W . (Al) 

Apart from the expression determining </>; ns t, we can 
retain most of the expressions from the first-order model 
and build a similar picture in the ng—r plane. We present 
this in Fig. [7J again for A = A' = 1. We see that the 
false vacuum dominated regime, which has ns > 1 and 
negligible r, lies entirely outside the WMAP5 allowed 
region, as does the main curve of the intermediate regime. 
Only once the trajectory heads towards the slow-roll limit 
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FIG. 6: Excluded regions in the (M, a) parameter space. The continuous (blue) line corresponds to ensuring model consistency; 
if the model lies above this region the phase transition will not take place. The long dashed (red) line corresponds to WMAP5 
constraints on the value of the scalar perturbations tilt and allows models to the left of the bound. The region between the 
short dashed (black) lines indicates models satisfying the maximum size of bubbles allowed by the level of anisotropy in the 
CMB. 



does it become compatible with observations. At M ~ and at M ~ 5 X 10 4 mpi inflation ends through slow roll 
9 x 10 _4 mpi the models cross the WMAP5 95% contour violation instead of a phase transition. 
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FIG. 7: Trajectory in parameter space ns — r describing second-order hybrid inflation models when M evolves from small 
through large values. The potential has one less coupling compared to the first-order case and all models are described by a 
single curve, as opposed to Fig. [3] Allowed models are those which approximate slow roll behaviour. False vacuum, blue tilted, 
models all lie outside the 95% C.L. 



